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Analysis of reinforced concrete deep beam is based on 
simplified approximate method due to the complexity of the 
exact analysis. The complexity is due to a number of 
parameters affecting its response. To evaluate some of this 
parameters, finite element study of the structural behavior of 
the reinforced self-compacting concrete deep beam was 
carried out using Abaqus finite element modeling tool. The 
model was validated against experimental data from _ the 
literature. The parametric effects of varied concrete 
compressive strength, vertical web reinforcement ratio and 
horizontal web reinforcement ratio on the beam were tested 
on eight (8) different specimens under four points loads. The 
results of the validation work showed good agreement with 
the experimental studies. The parametric study revealed that 
the concrete compressive strength most _ significantly 
influenced the specimens’ response with the average of 
41.1% and 49 % increment in the diagonal cracking and 
ultimate load respectively due to doubling of concrete 
compressive strength. Although the increase in_ horizontal 
web reinforcement ratio from 0.31 % to 0.63 % lead to 
average of 6.24 % increment on the diagonal cracking load, 
it does not influence the ultimate strength and the load- 
deflection response of the beams. Similar variation in 
vertical web reinforcement ratio leads to an average of 2.4 % 
and 15 % increment in cracking and _ ultimate load 
respectively with no appreciable effect on the load-deflection 
response. 
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1. Introduction 


Deep beams are members that are loaded on one face and supported on the other face such that 
strut-like compression element can develop between the loads and the supports provided that the 
clear span not exceed four times the overall depth or concentrated load exists within a distance of 
two times the overall depth of the beam from the face of the support [1]. Euro Code [2] on the 
hand does not account for the difference in the loading distribution. It defines a deep beam as a 
member whose span is less to or equal to 3 times the overall section depth. Due to the small span 
to depth ratio, deep beam supports heavy loads with little or no deflection, hence it is utilized in 
many structural applications requiring heavy load transfer such as bridges, diaphragms, water 
tanks, foundations, bunkers, girders used in multi-story buildings to provide column offsets and 
floor slabs under horizontal loads [3]. Typically, deep beams have a narrow width and contain a 
congested arrangement of reinforcements making it difficult for convectional vibrated concrete 
(VC) to flow easily through its web and adequately fill the bottom part. This often results in 
many problems in concrete such as voids, segregation, weak bond with reinforcement bars and 
holes in its surface [4]. It is on this ground that some studies have recommended the use of self- 


compacting concrete (SCC) for casting structural members such as deep beam [4]. 


SCC is an innovative concrete that is able to flow under its weight, completely fill the formwork 
and achieve full compaction, even in the presence of congested reinforcement without 
segregation [5]. It has gained preference over vibrated concrete due to a number of factors, 
including faster construction, reduction in site workforce, better surface finishes, easier placing, 
improved durability, greater freedom in the design of thinner concrete sections, reduced noise 
levels, the absence of vibration and safer working environment [6]. Comprehensive studies have 


been presented in the literature on the compositions and mechanical properties of SCC [5,7]. 


The difference in some of the mechanical properties and mix compositions of the vibrated 
concrete (VC) and the self-compacting concrete (SCC) demands an investigation of the structural 
behavior of members constructed using this type of concrete. Of major concern to researchers is 
the shear strength of deep beams. It is argued that the aggregate interlock mechanism, which is 
partly responsible for the shear strength of the structure, is likely to be weaker in SCC since it 
contains fewer and smaller sizes of coarse aggregate compare to vibrated concrete. Choi et al. [8] 
experimented with four deep beam specimens (two made with SCC and two with VC) to 
evaluate the influence of the concrete types on the shear strength of deep beams. The study 
reported that the SCC specimen having normal shear reinforcement distribution showed slightly 
higher performance than the corresponding VC specimen, while the SCC specimen having 


congested shear reinforcement condition showed a similar load carrying capacity to the 
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corresponding VC specimen. The work of Itterbeeck et al. [9] on shear capacity of SCC beams 
also reported improved shear strength over VC. On the contrary, Biolzi et al. (2014), as cited in 
[10], concluded that the ultimate shear strength of SCC beams is lower than that of VC beams. 


The structural design of deep beams is another area of great concern. The complex stress 
distribution in the deep beam cannot be resolved with the technical bending theory since the 
assumption of linear strain distribution cease to be valid and hence, the section of the beam 
belongs to disturbed-region (D-region) [11]. Previously, D-region has been designed based on the 
rule of thumb or purely experimental model [12]. The most recent approach to the design of D- 
region uses the strut and tie model (STM). It is a simplified method in which the continuous 
reinforced concrete domain is transformed to an equivalent truss as shown in Figure 1. The 
compression and tension zones are converted into equivalent struts and ties connected at the 


nodes to form a truss which can be easily resolved using basic mechanics [1]. 


p 


prismatic 
strut 


| Tie 


Fig. 1. Strut and tie model of deep beam [1] 


Despite the simplicity of the STM approach, there is no clear guide to define the geometry of the 
strut and tie model [13]. The major complexity is about transforming the continuous structural 
domain to optimal strut and tie model (STM). As reported by Liang et al. [13], Ritter (1899) 
found that a reinforced concrete beam after cracking due to diagonal tensile stresses could be 
modeled as a parallel chord truss with compressive diagonals inclined at 45° with respect to the 
longitudinal axis of the beam. Schlaich et al. [12] also recommended elastic stress method and 
load path method. 
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The reported techniques for developing STM involve trial and error iterative process which 
depend on the intuition and experience of the designer. Previous knowledge of stress distribution 
within the design domain is also necessary to reduce the number of iterative processes required 
to achieve adequate truss model. These lead to analysis result that is not unique but depends on 
the truss configurations. Therefore, the alternative design procedure is required to achieve 
uniform and efficient save design of deep beam. To achieve this, the proper characterization of 


various factors influencing the strength properties of the deep beam is worthy of evaluation. 


In the past, experimental procedure is mostly preferred as it provides real-life structural behavior 
that can be guaranteed of high accuracy. However, the need to execute the experimental study in 
line with standard specifications often make the procedure expensive and time-consuming. 
Besides, there are some other real-life situations that do not demand experimental work. In 
practice, there are many instances that performance assessment of structures incorporating deep 
beam is required due to changes in design codes, specifications, the occurrence of damage or 
degradation under service condition. Another case is when a constructed structure faces new 
conditions and needs to be retrofitted or repaired. Hence, the recent trends favored the use of 
numerical methods to evaluate the structural behavior of members. One of the most widely 
utilized approaches is the finite element methods (FEM). 


Ever since the introduction of FEM in Civil Engineering, it has been an effective tool for 
modeling and simulation of both linear and nonlinear behavior of reinforced concrete structures. 
Recent advancement in computational mechanics has also lead to the development of numerous 
material models that are compatible with finite element analysis method. In the present study, 
finite element model of the reinforced self-compacting concrete deep beam was developed using 
Abaqus finite element modeling tool. The study aimed to investigate the structural response of 
Reinforced Self-Compacting Concrete (RSCC) deep beam. The specific objectives of the study 
were (i) to validate the capability of the finite element (FE) model to capture the behavior of 
RSCC deep beam up till failure and (ii) to evaluate the parametric effect of some of the beam 
parameters on the structural behavior of the beam. The parameters studied are a concrete 
compressive strength, vertical web reinforcement and horizontal web reinforcement of the 


beams. 


2. Research program 


This study was conducted in two phases. The first phase of the work involved the validation of 
the FE model results versus the experimental values from the literature. The second phase 
focused on the parametric study to evaluate the influence of various design parameter on the 


load-deflection response, failure load, as well as the failure mode of the deep beam. 
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2.1. Model validation 


In order to demonstrate the capability of the FE model to simulate the response of reinforced 
self-compacting concrete (RSCC) deep beam under load up till failure, data on two (2) RSCC 
deep beam specimens previously tested by Choi et al. [8] were obtained from the literature. The 
beams have the similar geometrical properties but different shear reinforcement distributions; 
were tested to failure under 4-point loadings. The goal of the tests was to study the behavior and 
performance RSCC deep beam. The detailed configuration of the experimental set-up is shown 
in Figure 2 while the reported results of the mechanical property tests on reinforcing steels and 


concrete were presented in Table 1 and Table 2 respectively. 


414 
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i700 180 
(a) SCC-50 specimen 


360 


<I 
ra 


4700 
(b) SCC-100 specimen 


Fig. 2. Details and Dimensions of deep of beams (All dimensions in mm) [8]. 


Table 1 
Mechanical properties for the reinforcement steel [8]. 
Bar Diameter Yield Stress Ultimate Stress Elongation 


(mm) (N/mm”) (N/mm’) (%) 
25 334 512 33 
22 347 542 25 


6 447 743 21 
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Table 2 
Mechanical properties for the self-compacting concrete [8]. 
Compressive strength, f, Tensile Strength, f, Modulus of Elasticity, E 


(N/mm/”) (N/mm’) (N/mm/’) 


Mix Code 


HSCC py a | 4.2* 37677.76* 


The tensile strength, f; and the modulus of elasticity, E for the concrete were not reported in the 
study. The parameters were however obtained using the mathematical models in equations (1) 
and (2) respectively as given in [2]. 


fe = 0.3 * f,/3 (1) 
E = 22[(f- + 8)/10]°° * 103 (2) 


2.2. Parametric study 


The primary goal of the numerical simulation is to create a platform for parametric studies. In 
this work, the validated FE model was employed to evaluate the influence of various beam 
design parameters on the behavior of RSCC deep beam. Eight deep beams of the same 
geometrical dimension but different concrete compressive strengths, vertical and horizontal web 
reinforcement distributions as shown in Table 3 were modeled and analyzed. The beams’ cross- 
section dimension is 180 by 360 mm, and the effective length is 1300 mm. The reported material 


properties in Table 1 and 2 were adopted for the parametric studies. 


Table 3 
Parameters of the modeled deep beam specimens. 
Total ees Horizontal Pencil 
Beam L . i Section id Tension Compression — Vertical Shear , h xe Compressive 
Code sack Dimension * Reinforcement Reinforcement Reinforcement fe ies strength 
(mm) Reinforcement 2 
(mm) (N/mm‘) 
Bl 1700 180x360 1.4 2625 +2025 2022 6/100 06 / 100 52.10 
B2 1700 180x360 1.4 2625 +2025 2022 6/50 06 / 100 52.10 
B3 1700 180x360 1.4 2625 +2025 2022 6/100 06 / 100 26.05 
B4 1700 180x360 1.4 2625 +2025 2022 6/50 06 / 100 26.05 
B5 1700 180x360 1.4 2625+ 20625 2022 6/100 06 / 50 52.10 
B6 1700 180x360 1.4 20625 +2025 2022 6/50 06 /50 52.10 
B7 1700 180x360 1.4 2625 +2025 2022 6/100 06 / 50 26.05 
B8 1700 180x360 1.4 20625 +2025 2022 6/50 06 / 50 26.05 


42 M. A. Akinpelu, A. A. Adedeji/ Journal of Soft Computing in Civil Engineering 2-1 (2018) 36-61 


3. Finite element study 


The general purpose FE software, Abaqus [14] was used to develop FE models to mimic the 
structural response of the studied deep beams. The results of the simulation were compared to 
experimental results in term of load-deflection responses, diagonal cracking and ultimate loads 
as well as failure patterns of the beams. 


3.1. Concrete constitutive model 


Abaqus software offers three different constitutive models to simulate the nonlinear response of 
concrete namely (1) Smeared crack model; (2) Discrete crack model; and (3) Damage plasticity 
model [14,15]. The concrete damage plasticity model (CDPM) is the most preferred for this 
work due to the excellent performance recorded in similar studies [16,17]. The CDPM as 
implemented in Abaqus software is based on the models proposed by Lubliner et al. (1989) and 
Lee and Fenves (1998) as cited in [15]. The model combines the isotropic damaged elasticity 
model with isotropic tensile and compressive plasticity to represent the inelastic behavior of 
concrete. It assumes that the main two failure mechanisms are tensile cracking and compressive 
crushing of the concrete material. The evolution of the yielding (or failure) is controlled by 
tensile and compressive equivalent plastic strains under tension and compression loading 
respectively. Figure 3 illustrates the uniaxial tensile and compressive behavior of concrete as 
used in the CDPM. As depicted in the figure, the elastic stiffness of the material in the unloaded 
state is damaged or degraded. The evolution of damage in tension and compression is controlled 
by the introduction of tensile and compressive damage variables, d, and d. which takes on values 
between 0 and | for initially undamaged material and completely damaged material respectively. 
Further explanation of the CDPM model can be found in [14,15]. 


5. (b) Compression 


(a) Tension 


Fig. 3. Post-failure stress-displacement curve [14]. 
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3.1.1. Stress-strain relationship for concrete in compression 

The stress-strain data for concrete in compression and tension is crucial for analysis utilizing 
CDPM. However, this was not reported as part of the experimental results, only the ultimate 
compressive strength of the concrete at 28 days was recorded as given in Table 2. The curves 
were created using mathematical models found in the literature. The compressive stress behavior 
of concrete was generated using Eurocode 2 [2] model for non-linear structural analysis of 


concrete as shown in Figure 4. 


stress 


fom 


&. E, Strain 


Fig. 4. Stress-strain model for concrete in compression [2]. 


The relationship between the compressive stress and shortening strain for short-term uniaxial 
loading is as expressed in Equation (3). 


Oc _ _kn-1? 3) 
fom 1+(k-2)n 
The parameters in Equation (3) are expressed as follows: 
_ i 4 
1 (4) 


= O7ie (100 = 28 (5) 
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€, = 0.0035 (f; < 50N/mm?) (6) 
€y = (2.8 + 27[(98 — fem)/100] 4/100 (f > 50N/mm?) (7) 
k = 11E x (€y/ fem) (8) 
fom =f, + 8 (N/mm*) (9) 
Where: 


f, is the characteristic cylindrical compressive strenght of concrete at 28 days 
E, 1s the strain at peak stress 
E,, is the ultimate strain 


3.1.2. Tensile stress behaviour of concrete 

The nonlinear behavior of concrete in tension is characterized by a stress-crack displacement 
response instead of a stress-strain relationship due to its brittle behavior. This approach is 
preferred due to the mesh sensitivity problem associated with a stress-strain model for concrete 
in tension when there is no reinforcement in the significant region of the model [14]. The stress- 
crack displacement relationship can be model as linear, bilinear or exponential tension softening 
response as shown in Figure 5. For this study, the exponential stiffening curve according to 
Cornelissen et al., 1986, as cited in [18], was adopted. 


exponential 


----= linear 


— — bilinear 


Tensile Stress (MPa) 


0.8G,/f,  2G;/; 3.6G;/f, 5.14G,/F, 
W (mm) 


Fig. 5. Stress-crack width for concrete in tension [18]. 
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The model governing equation is as stated in Equations (10) and (11). 
o/fe= (1+ (22) Jem — (w/wJ(1 + eee (10) 
We = 5.145 (11) 


Cl and C2 are materials constants given as 3 and 6.93 for normal concrete, Wc is the critical 


crack opening width. 


Gr is the fracture energy per unit area required to create, propagate and fully break a Mode 1 


crack in concrete as given in Equation (12). 
Gp = Gro(fem/10)°’ (12) 


Gr is a coefficient which depends on the maximum size of aggregate. The coefficient can be 


obtained from [19]. The value used in this study is 0.08 J/m’. 


3.1.3. Other CDPM parameters 

(a) Dilation Angle, ~ accounts for the increase in plastic volume in concrete beyond the critical 
stress value. Different values of dilation angle for concrete are found in the literature. Lubliner et 
al. (1989) as reported in [15] gives Ww as 30 for concrete. However, recent studies have 
successfully used higher values of y . For instance, Sumer and Aktas [20] used 37° while Demir 
et al. [21] adopted 56°. Therefore, the appropriate value of w for this study was determined 


through sensitivity analysis. The values considered are 36°, 46° and 56°. 


(b) Viscosity Parameter, 1 is required for viscoplastic regularisation of the model. The value 
used for the study is 0.0001. The value is sufficiently small to enhance the rate of convergence 


without compromising the model results [22]. 


(c) Flow Potential Eccentricity, € defines the rate at which the surface approaches its asymptote 
(the flow potential tends to a straight line as the eccentricity tends to zero). The parameter is 


calculated as the ratio of uniaxial tensile strength to uniaxial compressive strength [23]. 


(d) The ratio of initial equibiaxial compressive yield stress to initial uniaxial compressive yield 
bo 


stress, =n was taken as 1.16 [15]. The ratio determines the point at which the material fails 
cO 


under biaxial compression. 
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(e) K Parameter is the ratio of the second stress invariant on the tensile meridian to that on the 
compressive meridian at initial yield for any given value of the pressure invariant such that the 


maximum principal stress is negative. The default value of 0.667 was adopted [15]. 


3.2. Steel reinforcement constitutive model 


The steel material is assumed to be linearly elastic. The elastic modulus of 210 N/mm? and 
poison ratio of 0.3 were adopted for the material. The two elastic constants define the isotropic 
elastic properties of the steel material. The experimental value of yield strengths presented in 
Table 1 with zero (0) plastic strain was used to define the plasticity property of the material. 
These values defined the steel materials that were applied to the truss elements that modeled the 


steel reinforcements in the FE model. 


3.3. Finite element model 


All the studied beam specimens were model in Abaqus. The concrete was modeled with an 8- 
node linear solid element with reduced integration to prevent shear locking effect associated with 
fully integrated solid elements as noted in [14]. The reinforcement steel was modeled with 2- 
node 3D truss element since it supports only axial loads. The bonding of reinforcement with the 
concrete was ensured with an embedded constraint with the concrete as the host region. The 
mesh sizes considered for the study are 25 mm, 50 mm and 75 mm in all directions. This was 
done to observe the influence of mesh density on the results of the model. For the boundary 
conditions, the nodes at the supports were restrained against vertical displacement while the 
external nodes along the symmetry planes were restrained along the longitudinal and transverse 
directions respectively to allow for Poisson’s effect. All the beams were loaded to failure by 
displacement control in the vertical direction from the datum point located above the beam and 
connected to the beam loading surface using equation constraint. Figure 6 shows a typical FE 
model for the beam section while Figure 7 shows the reinforcement model as embedded in the 
FE model. 
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Fig. 6. Typical FE model of the deep beam specimens (Specimen B2) 


Fig. 7. Typical FE model of the reinforcement for the beam specimens (Specimen B2). 


47 
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4. Results and discussion 


4.1. Model sensitivity analysis 

4.1.1. Dilation angle sensitivity analysis 

Different values of dilation angle were reported for concrete in literature hence, the need to 
calibrate the parameter to select the most appropriate for the studied beam specimens. Figure 8 
shows the results of the sensitivity tests. From the figure, the value of dilation angle that best 
modeled the experimental test is 56° and therefore used throughout the study. Table 4 gives the 
summary of other CDPM parameters finally adopted in the study. 


—s— FE Model (y=36°) 
—"— FE Model (w=46°) 


500 
a -—=— FE Model(w=56°) 
~ 400 —=— Experimental 
o 
oO 300 
as 
200 
100 
0 
0 1 2 a 4 5 6 z 8 
Deflection (mm) 
Fig. 8. Dilation angle sensitivity test. 
Table 4 
CDPM Model Parameters. 
Dilation Angle, p 30° 
Viscosity, 0.0001 
Eccentricity 0.1 
Stress ratio 1.16 


K 0.667 


M. A. Akinpelu, A. A. Adedeji/ Journal of Soft Computing in Civil Engineering 2-1 (2018) 36-61 49 


4.1.2. Mesh sensitivity analysis 

Generally, FE models are sensitive to mesh density. Figure 9 shows the effect of mesh density on 
the load-deflection response of the modeled beam specimen. The model with finer meshes 
captures the experimental result better than the model with coarse mesh although at the higher 
computational time. The best convergence result is obtained with a mesh size of 25 mm. This 


mesh size was used for the parametric study. 


1000 
900 
800 
700 
600 
500 
400 
300 
200 
100 


—— Experimental 
—s— Mesh 25 x 25 
—=— Mesh 50 x 50 
—=— Mesh 75 x 75 


Load (kN) 


0 1 2 3 4 5 6 vf 8 9 10 


Deflection (mm) 
Fig. 9. Mesh sensitivity test. 


4.2. Model validation 

4.2.1. Load deflection response 

The comparison of the load-deflection response of the FE model and the experimental results of 
the beam specimens are illustrated in Figure 10. From the figure, it can be inferred that the 
selected model is capable of capturing the experimentally observed loading trends and 
magnitudes for the entire loading regime. Although some deviations were noticed in the load- 
deflection response of the specimens, these could be attributed to several factors. The higher 


initial stiffness of the FE model results could be traced to the likely presence of micro-cracks in 
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the experimental specimen. These could occur through dry shrinkage of concrete and handling of 
the beams. It is also worthy to state that the accuracy of the materials parameters strongly 


influenced the observed response. 


Notwithstanding the above, the predicted failure loads and mid-span deflections compared well 
with the experimental values. As summarised in Table 5, the mean ratio of experimental to FE 


model results of the failure loads and deflections is 1.05 and 0.73 respectively. 


Load (KN) 


—*e— FE Model (SCC-50) 
—*— Experimental (SCC-50) 


0 1 2 3 4 5 6 7 8 9 
Defelection (mm) 


Load (KN) 


—e— FE Model (SCC-100) 
—*— Experimental (SCC-100) 


) 1 2 3 4 5 6 7 8 9g 
Defelection (mm) 


Fig. 10. FE model versus experimental load-deflection curve for the beams. 


M. A. Akinpelu, A. A. Adedeji/ Journal of Soft Computing in Civil Engineering 2-1 (2018) 36-61 51 


4.2.2. Diagonal cracking and reserved capacity 

At initial loading stages, elastic stress distribution exists in all the specimen leading to flexural 
cracks at the lower part of the beam between loading points in the mid-span region (Figure 11). 
As the loading increases, the first diagonal crack developed within the shear span; the elastic 
stress distribution was disrupted hence the load path changed, and forces flowed directly from 
the loaded points to the supports. This phenomenon is known as arch action in which 
compression struts connect the loaded points to the support and the tensile reinforcements act as 
ties linking the supports. The wide margin between the diagonal cracking strength and the 
ultimate strength could be linked to the arch action mechanism. The average ratio of the 
experimental to FE model reserved capacity of the beams, which is expressed as the ratio of the 
diagonal cracking load to the failure/ultimate load (Pc/Pu), is 1.01 as given in Table 5. 

Table 5 


Experimental and FE model diagonal cracking load (Pc), ultimate load (Pu), Reserved Capacity and 
maximum deflection (6 ) for RSCC deep beams. 


Experimental Experimental / FE Model 
Beam Reserve Reserve 
Cede Pc Pu Pc Pu 6 
Capacity oO (mm) Capacity 
(kN) (kN) (KN) (KN) (mm) 
(Pc/Pu) (Pc/Pu) 
SCC 
245 Yale: 0.34 2.10 101 1.01 1.00 0.68 
100 
SCC 50 324 796 0.34 2.60 1.10 1.09 1.01 0.78 
Average 1.06 = 1.05 1.01 0.73 


4.2.3. Failure mode 

The failures of reinforced concrete members are characterized by the occurrence and 
development of cracks, crushing of concrete as well as yielding of reinforcement. The crack 
patterns for the beams compare well with the experimental results throughout the loading regime. 
Figure 10 compares the cracks observe during the experiment with the cracks predicted in the 
simulations for SCC-50 specimen. The figure shows that at peak load, the beam failed through 
the diagonal crack developed within the shear span. Although there are some flexural cracks at 


the bottom of the mid-span, it does penetrate the compression strut. Therefore, the failure of the 
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beam is a shear failure. It needs to be emphasized that the concrete damaged plasticity model 
does not have the notion of cracks developing at the material integration point. However, it is 
possible to introduce the concept of an effective crack direction with the purpose of obtaining a 
graphical visualization of the cracking patterns in the concrete structure. The criteria adopted in 
this model are the assumptions that cracking initiates at points where the tensile equivalent 
plastic strain is greater than zero, and the maximum principal plastic strain is positive. Based on 
these criteria, the direction of the vector normal to the crack plane was assumed to be parallel to 
the direction of the maximum principal plastic strain. This direction was viewed in the 


visualization module of the post-processor of Abaqus CAE. 


PE, Max. Principal 
+5,.947e-02 
+5.451¢e-02 
+4,956e-02 
+4.460e-02 


+2.478e-02 
+1.982e-02 
+1.487e-02 
+9.911e-03 
+4.956e-03 
+0.000e+00 


| Fig. 11. Experimental and FE model crack pattern at failure for specimen SCC-50. 


4.3. Parametric study 


This section discussed the influence of the studied beam parameters on the diagonal cracking 
(Pc) and ultimate load capacity (Pu) of the deep beam using the FE model. The parameters 
considered are concrete compressive strength (f-), vertical web reinforcement ratio (ry) and 


horizontal web reinforcement ratio (r;). 


4.3.1. Effect of concrete compressive strength (f,) 
Table 6 illustrates the effect of f, on the cracking and ultimate load for all the modeled beam 
specimens. The percentage increase in ultimate load due to doubling the f, value ranges from 


46.34 % to 50.39 % (average increase of 49 %). The improvements slightly reduce with 
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increasing vertical and horizontal web reinforcement ratio. The improvement in diagonal 


cracking load due to doubling the fc value ranges from 39.75 % to 42.96 % (average increase of 


41.1 %). The percentage increase slightly drops with increasing horizontal web reinforcement. 


The ratios of cracking to ultimate load range from 0.43 to 0.52 for NSCC beams while the range 


is 0.41 to 0.49 for HSCC beams. These show that the strength ratio decreases with increasing fo. 


This is similar to the experimental result of Al-Khafaji et al. [4]. 


3 4 5 6 
Deflection (mm) 
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Fig. 12. Effect of f, on the load-deflection response. 
Figure 12 shows the influence of f, on the load-deflection response of the beam specimens. From 


the figure, the increase in (f,) value reduces the deflection for all the loading regime. This could 


be traced to the fact that increase in f, value results in higher modulus of elasticity which then 
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increases the flexural rigidity of the member and therefore lowers the deflection as previously 
observed by Al-Khafaji et al. [4]. 


Table 6 
Effect of f, on diagonal cracking and ultimate load. 
Normal Strength High Strength % Variation due to increase in 
% (f= 26.0 N/mm’) (f, = 52.1 N/mm’) ft 
Pc Pu Pc/Pu- Pc Pu Pc/Pu Pc Pu 
r= 0-31 223 451.3 0.49 318.8 681.4 0.47 42.96 50.99 
031 963 229 526.85 0.43 324 783.87 0.41 41.48 48.78 
r= 0-31 239 459.26 0.52 334 686.53 0.49 39.75 49.49 
0.63 9.63 244 536.65 0.45 345 785.34 0.44 41.39 46.34 


4.3.2. Effect of vertical web reinforcement ratio (ry) 

The effects of r, on the cracking and the ultimate loads is shown in Table 7. Increase in ry 
generally has a positive influence on the cracking and the ultimate loads of the modeled beam 
specimens. The percentage increase in ultimate load due to a doubling of r, ranges from 14.39 to 
16.85 % (average increase of 15 %). The percentage increment slightly reduces with increasing 
value of compressive strength (f,). The implication is that the effect of the increase in f, 
outweighs the influence of the increase in ry. This agrees with the experimental work of 


Mohammadhassani et al. [24]. 


The percentage increment in ultimate load increases with increasing r, for normal strength beam 
while the case is reversed for high strength beams. This further buttress the contribution of high 
strength concrete to the ultimate strength. Therefore, exceeding the minimum value of horizontal 


web reinforcement for high strength beams is less beneficial to the ultimate strength. 


The percentage increase in the diagonal cracking load varies between 1.63 % to 3.29 % (average 
increase of 2.4 %). The effect of r, on the percentage increment due to r, depends on the 
compressive strength of concrete. The percentage improvement slightly drops with increasing rp 


for normal strength beam while it increases significantly with increasing r, for high strength 
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beam. This means that higher horizontal web reinforcement benefits more in crack control for 
high strength beams which is contrary to the conclusion of Mohammadhassani et al. [24] that the 


web reinforcements do no influence the diagonal cracking load. 


The ratio of cracking load to ultimate loads also decreases with increasing r,, this implies that 


increasing r, contribute more to the ultimate load than the cracking load. 


The plots in Figure 13 illustrate that the increase in vertical web reinforcement (ry) does not 
influence the load-deflection response of the beam. Therefore, each pair of beams that only differ 
in ry have a convergent load-deflection response. This is because ry does not contribute to the 


flexural rigidity of the beam as previously observed by Al-khafaji et al. [4]. 
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Fig. 13. Effect of r, on the load-deflection response. 


56 M. A. Akinpelu, A. A. Adedeji/ Journal of Soft Computing in Civil Engineering 2-1 (2018) 36-61 


Table 7 
Effect of r, on diagonal cracking and ultimate load. 


% Variation due to 


ry%= 0.31 ry% = 0.63 ; : 
increase in ry 
tn % 
Pe Py P./P, Per Pu Po/Pu Pe Py 
0.31 223 451.3 049 229 526.85 0.43 2.69 16.74 
f,= 26.0 
N/mm” 
0.63 239 459.26 0.52 244 536.65 0.45 2.09 16.85 
0.31 318.8 681.4 0.47 324 783.87 0.41 1.63 15.04 
f= 52.1 
N/mm” 
0.63 334 686.53 0.49 345 785.34 0.44 3.29 14.39 


4.3.3. Effect of horizontal web reinforcement ratio (1) 


The parametric effects of r, on the cracking and the ultimate loads of the beam specimens is 
computed in Table 8. Unlike vertical web reinforcement (r,), increase in r, only contributes little 
or no influence to the ultimate loads of the modeled beam specimens. The percentage increase in 
ultimate load due to a doubling of r, ranges from 0.19 % to 1.86 % (average increase of 1.14 %). 
The percentage increment reduces with increasing value of compressive strength (fc). This effect 
is similar to that of vertical web reinforcement. This result further established the influence of f, 


on the ultimate strength of the deep beam. 


The percentage increment in ultimate load increases with increasing r, for normal strength beam 
while the case is reversed for high strength beams. This is also similar to the result of horizontal 
web reinforcement (r,) although the vertical web reinforcement contributes more to the ultimate 


strength when compared to horizontal web reinforcement. 


As seen in table 8, increase in horizontal web reinforcement ratio generally improves the 
diagonal cracking load of the beams. This is also contrary to the observation of 


Mohammadhassani [24]. The percentage increase in the diagonal cracking varies between 4.77 
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% to 7.17 % (average increase of 6.24 %). The percentage improvement decreases with 
increasing f, for all the beam specimen. Also, the percentage improvement slightly drops with 
increasing r, for normal strength beam while it increases significantly with increasing ry for high 
strength beam. This means that combination of high horizontal and vertical web reinforcement is 


more effective in crack control for high strength beams than normal strength beams. 


As seen in Figure 14, variation in horizontal web reinforcement (r,) does not influence the load- 
deflection response of the studied beams. Therefore, each pair of beams that only differ in rj 
have convergent load-deflection response. This is also because r, does not contribute to the 


flexural rigidity of the beams. 
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Fig. 14. Effect of r, on the load-deflection response. 
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Table 8 
Effect of r, on diagonal cracking and ultimate load. 


% Variation due to 


th %= 0.31 th %= 0.63 ; ; 
increase in rp 
r% 
Pe P, P./Py Per Pu Po/Pu Ps Px 
0.31 223 451.30 0.49 239 459.26 0.52 7.17 1.76 
f= 26.0 
N/mm” 
0.63 229 526.85 0.43 244 536.65 0.45 6.55 1.86 
0.31 318.8 681.4 0.47 334 686.53 0.49 4.77 0.75 
f= 52.1 
N/mm? 
0.63 324 783.87 0.41 345 785.34 0.44 6.48 0.19 


5. Conclusions 


Finite element analysis of reinforced self-compacting concrete deep beam was conducted using 
Abaqus modeling tool. The model was validated against experimental results from literature and 
sensitivity of the model to varied dilation angle was studied. The parametric effect of varied 
value of concrete compressive strength, vertical web reinforcement ratio and horizontal web 
reinforcement ratio on the diagonal cracking and ultimate load capacity as well as the load- 
deflection response of deep beam were also evaluated. Based on the results of this study, the 
followings are concluded: 

(1) The FE model is capable of simulating the entire loading response of a reinforced 
concrete deep beam with reasonable accuracy. The mean ratio of experimental to the 
predicted value of diagonal strength is 1.06; ultimate capacity is 1.05 and mid-span 
deflection is 0.73. 

(2) The compressive strength of concrete has the most significant influence on the 
deflection, the cracking and the ultimate load of the deep beam when compared with 


other parameters investigated. This observation conforms with the experimental study of 
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Al-Khafaji et al. [4] and Mohammadhassani el at. [24] on self-compacting concrete deep 
beam. 

(3) The average percentage increase in diagonal cracking and ultimate load due to doubling 
of concrete compressive strength is 41.1% and 49 % respectively for the studied beams. 

(4) The vertical web reinforcement and the horizontal web reinforcement does not influence 
the load-deflection response of RSCC deep beam. 

(5) The vertical web reinforcement ratio (ry) of deep beam contributes more to the ultimate 
strength than to the diagonal cracking strength. The average percentage increase in the 
diagonal cracking and ultimate load due to increase in ry from 0.31 % to 0.63 % is 2.4 % 
and 15 % respectively. 

(6) The horizontal web reinforcement ratio does not substantially influence the ultimate 
strength of the RSCC deep beam although it has more influence on the diagonal cracking 
strength when compared with vertical web reinforcement ratio. The average increment in 
diagonal cracking strength is 6.24 % due to doubling of r, ratio while it is 2.4 % due to 
doubling of ry. Therefore, r, contributes more to the diagonal cracking strength while r, 


contributes more to the ultimate strength. 
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